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, Abstract 

j^Q, A description of the short time behavior of solutions of the Allen-Cahn equation with a 

■ smoothened additive noise is presented. The key result is that in the sharp interface limit 

' solutions move according to motion by mean curvature with an additional stochastic forcing. 

This extends a similar result of Funaki |9] in spatial dimension n = 2 to arbitrary dimensions. 

in 

^7 1 Introduction and main result 

^"J \ l. Setting and main result: For a small parameter e > consider the following stochastic Allen- 

Cahn equation in an open domain D in R™ for some n > 2: 

S3 ' d 

— u £ (x, t) = Au £ (x, t) + e~ 2 f(u(x, t)) + £ _1 £ e (f) (x, t) e D x [0, oo) 

ot 

—u £ (x,t) = Q xedD (1.1) 

ov 



o 



u e (a:,0) = Mq( x ) x £ D. 



> 

o 

00 . 

ly-j , Here f(u) — —F'(u) is the negative derivative of a symmetric double-well potential. For fixing 

ideas, assume that F(u) = — and f(u) = u — u 3 . In particular F has two global minima at 
±1 and solutions of the dynamical system x = f(x), that start outside of zero, converge to one of 
these minima. The expression t; £ (t) denotes a noise term defined on a probability space (fi, P). 
The noise £ £ (t) is constant in space and smooth in time. For e J. the correlation length goes to 
zero at a precise rate and J* £, £ (s)ds converges to a Brownian motion pathwisely. The details of 
the construction and further properties can be found below. 

We study the short time evolution of developed surfaces for Ql.lj) . More precisely let So the 
boundary of a set Uq be compactly embedded in D of class C 2 - a for some a > 0. Assume that 
the initial configuration u £ (x,t) is close to —1 on Uq and close to +1 on D \ Uq with a transition 
layer of order 0{e). We show that for short times there exist two phases and the evolution of the 
phase boundary follows two influences - the tendency to minimize the boundary and a stochastic 
effect. The main result is: 

Theorem 1.1. Consider the problem (COP with the noise term £, £ {t) as constructed below. In 
particular suppose that the approximation rate 7 verifies 7 < | . Then for any compactly embedded 
hypersurface So = dU^ of class C 2 ' a there exist initial conditions u £ , a positive stopping time r 
and randomly evolving closed hypersurfaces (S(t))o<t< T such that the following hold: 

(i) The surfaces (S t )o<t< T evolve according to stochastically perturbed motion by mean curva- 
ture, e.g. the normal velocity V at each point is given by 

V = (n-l)K-coW(t). 
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(ii) sup 0<t<T \\u £ (x,t) — XE t |U 2 (i5) almost surely as e goes to zero. 
Here k denotes the mean curvature of the surface at a given point. The constant Co is given by 

V2 



Co = 



The function xs t ls a step function taking the value —1 in the interior and +1 on the exterior. 
The precise meaning of the geometric evolution will be given in the next section. 

The noise scaling £ _1 £ £ (i) can be interpreted as follows: Consider the stochastic equation 

|=A, + /(,)+ E f(t). (1.2) 

Equation (|l.ip can be obtained from this equation by diffusive scaling: u(x,t) — v(s~ 1 x,e~ 2 t). 
The intuition is that in (|1.2p surfaces should move with velocity V — (n — l)n + c(e£ £ (i)). Here c 
is the speed of a travelling wave solution corresponding to a perturbation of the potential through 
s£, £ (t). Then after rescaling one obtains as normal velocity V = n + e~ 2 x e c(e£ e (t)) such that 
the random term becomes a quantity of order 0(1), The significant observation is that the noise 
term does not rescale. Actually this observation is characteristic for our result. Even in the limit 
the Brownian motion can be considered pathwise and there is nowhere any need to work with 
stochastic integrals. 

2. The white noise approximation: Let (W(t),t > 0) be a Brownian motion defined on a 
probability space (f2, T, P). For technical reasons extend the definition of (W(t), t > 0) to negative 
times by considering an independent Brownian motion (W(t),t > 0) and setting W(t) — W(—t) 
for t < 0. Then (W(t),t € R) is a gaussian process with independent stationary increments and 
a distinguished point W(0) — a.s. Let p be a mollifying kernel i.e. p: EL — > EL+ is smooth and 
symmetric with p(x) = outside of [—1, 1] and J p(x)dx — 1. For 7 > set p s {x) = e ' p(-j^). 
Then the approximated Brownian motion W £ (t) is defined as usual as 

/oo 
p e {t - s)W(s)ds. 
-00 

Note that it is only here that the Brownian motion at negative times is needed. So actually 
only negative times in (— e 7 ,0] will play a role. The parameter 7 determines how quickly the 
approximations converge to the true integrated white noise. We will always assume 

2 

in order to have the needed pathwise bounds on the white noise approximations. 

Proposition 1.2. Let^ E {t) = W e (t) denote the derivative ofW e . Then the following properties 
hold: 

(i) £, e (t) is a stationary centered gaussian process with E[£ e (t) 2 ] = e~ 1 \p\ 2 L2 - 
(ii) The correlation length of^ e (t) is 2e 7 i.e. if\s — t\ > 2e 7 then^ e {t) and^ £ {s) are independent. 
(Hi) If j < 7 for every positive time T there exists a non-random constant C such that 



P 

In particular for 7 < I for e small enough 



3s a s.th. Ve < eq sup \^(t)\ < Ce~? 

Q<t<T 



= 1. 



£ e (*)<Ce~ 3 . (1.3) 
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Proof. One can write 

where the first equality follows from differentiating under the integral and the second from stochas- 
tic integration by parts. Then properties (i) and (ii) follow from standart properties of the stochas- 
tic integral. To see (iii) write 



\e(t)\ 



< 



t+e 

t—en 

2 7 



t+e 



, ft — s 



t+e- 1 , I i 



W(s)ds 
W(t)ds 



r-27 



t- S 



(W{t) - W(s))ds 



The first term vanishes due to J t _ £ ~, p' (^pr) ds — 0. One obtains 



if wi < 



-2- 



t+e< 



t- S 



(W(t) - W(s))ds 



< 



1 - 2 ~*2e~<\\p'\\ OC! osc se[t _ £ - lyt+£ - l] W(s) . 

The oscillation is defined as osc se[t _ eT it+s7 ]W(s) := sup se[t _ £T t+e7 ] - inf s£ [t_ e -y, t+E -Y] W(s). 

Now one can apply Levy's well known result on the modulus of continuity of Brownian paths 
(See e.g. [10] Theorem 9.25 on page 114): 



P 



limsup— max \W(t) — W(s)\ = 1 
A -^o g(6)o<s<t<T< 

t-s<6 



where the modulus of continuity is given by g(S) — y251og(|). In particular there exists almost 

surely a (random!) e such that for e < £ we have sup te [ .T] osc s£[t-e~< ,t+e~<]W(s) < (2s 7 ) 3 - ^ . 
This gives the desired estimate 



|£ s (t)| < £-^211^1^(2^)5-^ < ce^l\ 



We will need a similar bound on the derivatives of £ £ 



□ 



Proposition 1.3. Consider the process £ e (i). Then ifj < 7 for every positive time T there exists 
a constant C such that 



P 



3e Ve<e sup \?(t)\ < Ce~^ 

0<t<T 



1. 



In particular for 7 < § and e small enough 

If Wl < cs- 1 

Proof. The proof is similar to the one above: 

't-s 



(1.4) 



< 



t-el 
t+e'< 



e-^p" 



e-^p" 



t - S 



W{s)ds 



{W{t) - W(s))ds 



< 2e-^\\p"\\ ooOS c se[t _ e i tt+ ^ ] W(s) : 
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Then one applies Levy's modulus of continuity again to see that almost surely for e < eo{u>) one 
has osc sS [ a uW(s) < (2e 7 ) 5 and obtains the desired result: 

t(t)\ < 2e-^\\p"\\ 00 (^-^ h - = Ce%. 

□ 

3. Motivation and related works: Solutions of the Allen-Cahn equation 

du . 1 „. . 
Tt = An +-/(„) 

evolve according to the 1? gradient flow of the real Ginzburg-Landau energy functional: 



n e (u) = 




There are two different effects. The reaction term e~ 2 f(u) pushes solutions to the two minima 
±1 and the diffusion term Am tends to smoothen the solution. For small e there will be two 
phases, corresponding to regions where the solution is close to ±1. The width of the transition 
layer between those two phases is of the order 0(e). Then the evolution gradually shrinks the 
transition layer. 

This behavior is the motivation to consider the Allen-Cahn equation as a simple model of a 
two phase system which is driven by the surface energy without conservation of mass. Allen and 
Cahn p] introduced it to model the interface motion between different cristaline structures in 
alloys. In the deterministic setting there were major advances in connection with the improved 
understanding of the theory of geometric flows of surfaces as initiated for example by [fH [2] in 
the early nineties. In particular in [5] it was shown that in the limit e J. solutions only attain 
the values ±1 and the phase boundary evolves according to motion by mean curvature. The key 
difficulty here is to find a description of the geometric evolution which is global in time. A similar 
result for short times was established in [T3] . 

Stochastic perturbations of this effect have also been considered. From a modelling point of 
view an additional noise term can account for inaccuracies of the simplified model or as effects of 
thermal perturbations. From a mathematical point of view it is a very interesting and challenging 
question to study stochastically perturbed evolutions of surfaces and the Allen-Cahn setup is one 
possible point of view. In [8] Funaki considered the case of the Allen-Cahn equation in one space 
dimension with a space-time white noise. He showed that in the limit e J, on the right time-scale 
solutions only attain values ±1 and the boundary point essentially performs a Brownian motion. 
In [9] he studies the two dimensional case with a smoothened noise and shows that for short 
times solutions evolve according to a stochastically perturbed motion by mean curvature. His 
analysis relies on a comparison theorem which requires the noise to be smooth and a very subtle 
analysis of a quasi-linear stochastic PDE which describes the boundary evolution. On the level 
of stochastic surface evolution there were advances by Yip [15] and Dirr, Luckhaus and Novaga 
[4] but a fully satisfactory description is not yet available. Some results based on a stochastic 
version of the concept of viscosity solutions were announced in [1.2]. Recently the model has 
enjoyed an increasing interest in the numerical analysis community. For example in numerical 
approximations of the one-dimensional equation are studied. Numerical analysis of this equation 
is challenging because all the interesting dynamics happen on a very thin layer which requires to 
develop adaptive methods which work in the stochastic setting. 

Our result is a generalization of Funaki's result to arbitrary dimension. We use the same 
comparison technique to study the equation. Therefore we also need to assume a smoothened 
noise with correlation length going to zero as e goes to zero. The description of the surface and 
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the convergence result is based on 0] and fully avoids Funaki's result of weak convergence. In fact 
this is also a strictly pathwise result so that all results hold almost surely. 

4- Structure of the paper: In Section 2 the technique of [4] to describe motion by mean cur- 
vature is briefly reviewed and the main results are stated. In Section 3 the results about the 
geometric flow are used to proof the behavior of the Allen-Cahn equation. 



Acknowledgement: The author expresses his sincere gratitude to Tadahisa Funaki for the 
great hospitality he received at the University of Tokyo. He also thanks the referee for careful 
reading and various suggestions. 

2 Stochastic motion by mean curvature 

This section reviews the description of a stochastically perturbed motion by mean curvature given 
in [4]. A short time existence result for surfaces moving with normal velocity dV^ = (n — l)/tdt + 
cdW(t), where k denotes the mean curvature, and a pathwise stability result under approximations 
of the integrated noise are given. 

Motivated by [7] consider the following system 

dd(x, t) = g(D 2 d(x, t),d(x, t))dt + dW(t) (x, t) G O x (0, T) 

|Vd| 2 = l (x, t) e dO x (0, T) (2.1) 

d(x,0)=d o (x) xeO, 

on some open bounded domain O. Here D 2 d denotes the Hessian of d and g(A, q) = tr(A(i— qA)^ 1 ) 
for a symmetric matrix A and g£E. The initial condition do is supposed to be of class C 2a and 
to verify |Vd| = 1 in O. Furthermore it is assumed that Vc? is nowhere tangent to the boundary. 

In order to solve the above system consider q{x,t) = d(x,t) — W(t). Then q solves the system 

dq(x, t) = g(D 2 q(x, t), q{x, t) + W{t))dt (x, t) £ O x (0, T) 

|Vg| 2 = l (x, t) £ dO x (0, T) (2.2) 

q(x, 0) = d (x) x E O. 

Due to maximal regularity of the linearized system ([13]) and a fix point argument the following 
results are obtained: 

Theorem 2.1. ([1] Section 4) Let t 1— > W(t) be a-Holder continuous for some a 6 (0, 1). Then 
there exists a time T depending only on the C a l 2 -norm of W and the C 2 ' a -norm of do such that 
on O x [0,T] system $2.2\) and therefore also (2.1]) admit a unique solution of class C 1+a / 2 ' 2+a . 
Moreover if 1 1— > W{t) is another function of class C a and q is the solution \2. 2]) with W replaced 
by Wwith interval of existence [0,T] on has 

sup \\q(t,-)-q(t,-)\\c*,« <C\\W-W\\ Ca/H[0<raill{iiT}] y (2.3) 

tG[0,min{f ,T}] 

Now let £0 = dUg be as above. In particular So is assumed to be of class C 2,a . Define the 
signed distance function dg and the indicator xs as 



d (x) 

and 



Xs (x) 



— dist(a:, So) for x e Uq 

dist(cc, So) for x S D \ Uq 

for x £ Uq 
for x £ D \ U . 
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There exists an open environment O of So such that on O the function d (x) is of class C 2 ' a 
and Vd is nowhere tangent to dO. Furthermore on O it holds |Vdo| = 1- Then for a given 
stochastic noise W(t) consider the pathwise solution d(x, t) of l|2.ip with initial condition do on 
[0, T(lj)]. Define the evolving surfaces (£(t),0 <t< T(u)) as the zero level sets of d(x,t). Then 
the following holds: 

Theorem 2.2. ( g] Section 4) 

(i) For every t the function x i— » d(x,t) is the signed distance function o/S(t) on O. 

(ii) If X(0) in S(0). Then up to a stopping time there exists a solution X(t) to the stochastic 
differential equation 

dX(t) = (n - l)v[X(t),t)K(X(t), t)dt + v(X(t),t)dW(t), 

with X(t) G £(i) almost surely. 

Here v{x,t) denotes the exterior normal vector to £(i) for x £ £(i). The last observation 
justifies to say that the surfaces £(t) evolve according to stochastic motion by mean curvature. 
Note that we use the convention that K = — j J27=i Ki w ^ tn tne P rmc ip a l curvatures /tj such that 
the factor (n — 1) appears which is not present in [4]. 



3 Construction of sub- and supersolutions 

In this section the link between the boundary dynamic and the Allen-Cahn equation is established. 
For a related calculation see [9j E|- 

In order to construct sub- and supersolutions to (|l.ip consider the following modification of 
the reaction term: f(u, S) = f(u) + S. The implicit function theorem implies that there exists 
an interval [— ^o,^o] such that for 5 e [ — <5o, 5 ] there exist two solutions m±(5) of the equation 
f(u,6) — which are close to ±1 and that the mappings 8 i— * m±(8) are smooth. Consider the 
following auxiliary one dimensional problem 

d 0^ 

—u(x, t) = g^u(x, t) + f{u{x, i)) + 8 (3.1) 
u(±oo) = m±(8). 

A travelling wave solution to (|3.1|) is a solution u{x,t) = m(x — ct) with a fixed wavespeed c. 
Finding such a solution is equivalent to finding an appropriate waveshape m(x, 8) and wavespeed 
c(<5) such that 

m"(x) + c{8)m'{x) + {/(m(x)) + 8} = (3.2) 

m(±oo) = m±(8). 

The following properties hold: 

Lemma 3.1. ([3] Lemma 3.3) There exists a constant So such that for 8 e [— <5o,#o] problem 
$3.2\) admits a solution {mix, 8), c{8)) where m is increasing in x and this solution is unique up to 
translation. Furthermore m can be chosen smooth in 8. There exist constants A and (3 such that 
the following properties hold: 

(i) < d x m(x, 8) < A for all (x, 8) e R x [-S , do] ■ 

(ii) \d x m(±x,S)\ + \(d x ) 2 m(±x,8)\ + \m{±x,S)-m ± (S)\ < Ae~ 0x for all (x,S) e R+ x [-S ,6 ]. 
(Hi) The traveling wave velocity c(S) is smooth in [—80,80] and c(0) = 0. 
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Actually as pointed out in [9] dsc(0) = — c = — ^ — y= 

The idea of the construction is the following: We expect the surface to evolve according to two 
influences - the surface tension and the stochastic perturbation of the potential making one of the 
stable states more attractive. Close to the surface the solution should look like a travelling wave 
interface which is moving with velocity c(e£). This means that solution should behave like 

u(ar,t)ram(^^W(*) 



where d is the signed distance function of a surface moving with normal velocity V = (n — 1)k + 
e^ 1 c(e^ e ). The standard way of making this idea rigorous is to modify it in such a way that such 
an approximate solution is a true sub/supersolution and show that the difference between the two 
cases evolves on a slower time scale than the original dynamic. 

Fix some initial surface S as in Theorem ll.il As So is compactly embedded one can fix an 
N such that all the principle curvatures of So are bounded by N. As in Section [2] one can define 
a random evolution (S ±,£ (i),0 < t < T^ N ) evolving with normal velocity 

V = {n- 1)k + e^cie^it) ± e 13 ). 

Here the stopping time T^ N is denned as the largest time such that the evolution is well denned 
and such that on [0, T e N ] the principle curvatures remain bounded by N. The constant (3 can be 
chosen such that 1 < f3 < 2, The condition (3 > 1 ensures that in the original time scale the extra 
term does not have an effect and the condition (3 < 2 ensures that the effect is strong enough for 
the solution to remain a sub/supersolution. Furthermore assume (by shortening the time interval 
if necessary) that there exists an open set O such that for all t £ [0, T £ N ] the 77-neighborhood 
of S ± (t) is contained in O for some small 77. Then one can extend the signed distance functions 
d ± (x, t) to a smooth function on all of [0, T(w)] x D such that on U(t) \ O the function cP is 
smaller than —77 and on D \ (U(t) U O) it is larger than rj, such that iVc^l < 1 and such that d is 
constant close to 3D. 



Define 



±, x / d ± (x,t) ±e a e Clt , , B \ 



£ 

where a and c\ are constants that will be chosen below. One gets the following conclusion: 

Lemma 3.2. // one chooses a and c\ properly, there exists a (random) Eq > such that for all 
e < £0 and for < t < T N ' e 

U e '~{x,t) < M £ (x,t) < U e ' + (x,t), 

for every solution u £ (x,t) of (1.1)) with initial data verifying u e, ~(x,0) < u £ (x,0) < u e ' + (x,0). 

Proof. The conclusion will follow by a PDE-comparison principle. We only show the inequality 
involving u + the other one being similar. Let us calculate 

d t u E '+(x,t) = ^(d t d(x,t) + e a c ie Clt ^ + em s £(t) 
Au E (x, t) = ^Ad(x. t) + ^\Wd(x, t)\ 2 . 

Here m x denotes the partial derivative of m{x, t) with respect to x. Then rewrite the reaction 
term using p.2p : 

£- 2 {f{m) + ee e ) = e- 2 f-m" - m'c{e£, e + e") - 
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By properly arranging the terms one gets 

C(u+) := d t u e >+(x, t) - Au e (x, t) - £- 2 (/(m(x, <)) + e^(t)) = h + h + h + £ ~ 2 , 

where 

h = ^j- (d t d(x, t) + e a c ie Clt - Ad(x, t) + e- 1 ^ (t) + e* 3 )) 
h = emsiit) 

Here the first term accounts for the boundary motion. The statement that this term is small 
essentially means that the surface evolves with normal velocity V = (n — + e~ x c{e^ + e 13 ). 
The second term corresponds to the change of wave profile due to the change of noise. It is here 
that we need the pathwise bound (|1.4p on the derivative of £ e to control this term. The third term 
essentially vanishes because close to E(t) the function d coincides with d and therefore |Ve?| 2 = 1. 
Off the boundary the derivative m xx becomes exponentially small such that we also control this 
term. In the end this means that the correction term e@~ 2 dominates the dynamic. Let us make 
these considerations rigorous: 

By l|1.4p I2 < C for every e smaller than £q(ui). For d(x,t) < r\ Vd(x,t) = 1 and therefore I3 
vanishes for such x. For d(x,t) > r\ Lemma I3TT1 (ii) implies: 

2_A 

7 2 

To bound I\ consider points x close to £(t). For all other x the reasoning is as for I 3 . For x with 
dist(x, £ ±,e (i)) < 5^ the functions d{x,t) and d(x,t) coincide and one obtains 

d t d(x, t) = Ad(y, t) + e- 1 c(ee £ (t) + s ), 

where as before y is the unique point in £(t) such that d(x,t) = dist(x,y). Plugging this into I\ 
gives 

h = ^f {Ad e (y,t)~ Ad e (x,t)+e a Cl e Clt y 

Here one uses the fact that all the principle curvatures Ki{t,y) of the H ±,e {t) are bounded by N 
to obtain 



F(l"|Vd| 2 ) <M e -C/e^ _ 



\Ad E (y,t) - Ad e (x,t)\ 



n — 1 n—1 , ,\ 

^' n l-d(x,t)\ Kl (y,t)\ 
< 4N 2 d(x,t), 

because sup Ig [ ^i] 9 x j^ = 4. Plugging this in yields 

11,1 < m x {t^±^,ee±e^ 

Choosing c\ larger than A*" 2 and using sup,,, xm x < 00 one obtains |7i| < C. Thus altogether if e 
is small enough the term e' 3 " 2 will dominate everything else and one obtains 



£(«+) > 0. 
ition of d. S( 

desired result. The inequality for u- is shown in a similar manner. □ 



On the boundary = due to the definition of d. So a standart comparison principle gives the 
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To finish the proof of the main theorem one needs the following Lemma: 

Lemma 3.3. Fix any time interval [0,T]. Denote by W ± ' £ the random functions [0,T] 3 t 
^ f Q e~ 1 c(e^ e (s) ± e /3 )ds. Then c W ± ' E converges almost surely tot^ c W(t) in C°' Q ([0, T}) for 
every a < h. 

Proof. Consider only W + ' £ (t) the calculation for W~' s (t) being the same. Fix a < \ and a d with 
a < d < \. Then for P-almost every ui there exists a random constant C such that 

\W(s)-W(t)\ 



sup 

-l<s<t<T 



s-t 



< C. 



Assume that e is small enough to ensure s£ £ (t) + G [— 5 Q , S ]. (Recall that c is only defined on 
[—Sq, So] ) Using Taylor-formula and c(0) = one can write for every t: 

e- l c{eC{t) + eC) = c'(0)(£ e (i) + e^ 1 ) + l c "{a{t))e- x (ef (i) + e?)\ 

for some a{t) verifying \a{t)\ < \s£, e (t) +e /3 |. Therefore one can write 

||coT^-co^ + ' e ||oo < sup \coW(s)-c [ e^dtl+coTe^ 1 
selo.T] Jo 



£f„\\2 _L =-2/3-1 



+ T sup \c"(S)\{ sup e^(s)Y+s 
se[-s„j ] y se[o,T] 

Due to (|1.3p the last terms converge to zero almost surely. Therefore it remains to consider the 
first term. Due to W e (s) = £ e (s) one obtains: 

sup \coW(s)-co[ e(t)dt\ <cq sup \W(s)-W{sy\+c \W(Oy\ 



s6[0,T] 



s£[0,T] 



Co sup 

*6[0,T] 



(\V(s) - W{s-t)^ P e {t)dt +c j (w(0)-W(tj)p e (t)dt 



< 2c, 



0. 



Consider now the Holder-seminorm 
1 



sup 



0<s<t<T (t - s) a 

— sup 

0<s<t< 

< sup 



n > / T (t — S) 
1 



0<s<t<T {t — s) 
1 

SUp 7 r— 

0<s<t<T (t - S) a 



coW(t) - c W + - e (t) - coW(s) + coW+' s (s) 

c Q W{t) - coW(s) - J e^cie^is) + e /3 )ds 
\c W(t) - c W(s) - c W e (t) + coW £ (s)\ 



C £ 



0-1 



sup 

8e[-S„,5„ 



2 +£ 2/3-l 



du 



Again the second term converges to zero. For the first term one gets: 
1 



sup 



0<s<t<T (t - s) a 



\c Q W(t) - c W(s) - c W £ (t) + c W e (s)\ 



< c sup 

0<s<t<T (t — S) 



(W(t) - W(t -u)- W{s) + W(s - u))p £ (u)di 



< Co sup 

0<s<t<T 



(2(W(t) - W(s)))f 

< c (2C7)- (2C(2e' ( f) 1 ~ f . 
This shows the desired convergence. 



(W(t) - W(t - u) - W{s) + W(s -u)) s p e {u)du 



□ 
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Proof, (of Theorem Chose the initial configurations Uq such that u £ (x, 0) < — 

u e,+ (a;,0). Define the stopping time t{w) :— mf e T s N , Remark that r is almost surely posi- 
tive due to the boundedness of the || W^ ±,e || C c,/2 and the C 2,a convergence of d e,± to d. 

Then by Lemma [3.21 one has for all times < t < T{u>) that u e {x,t) < u £ (x,t) < 
So one gets: 

IK(-,t)-Xs (t )M)ll^ < \\u £ (- 7 t)-u e ' + {-,t)\\ L 2 {D) + \\u^+(-,t)- X ^, + (t)(-,t)\\ L 2 {D) 

+ ||Xs«.+ (*)(->*) - Xs(i)0>*)IU 2 (z>) 

< \\u £ '-{-,t) - U e >+(;t)\\ L 2- (D) + \\u e >+(;t) -X&.+ (*)(', t)\\ L 2 (D) + 

||XE«.+(t)(*)*) - XE(t)(-)*)IU 2 (r>) 

< \\u £ '-(-,t) - x ^--(t)(-,t)\\ L 2(D) +2\\u e -+(-,t) - x ^.+ (t)(;t)\\mD) + 

2 ||X£=.+(t)M) - X£(t)M)IU 2 (£>) + ||XS«. -(*)(■,*) - XE(t)(-,<)||i»(Z»)- 

The supremum in time of the first two terms converges to zero due to the definition of u e,± . 
Consider \\xv°,-(t)(') - Xs(t(-)IU 2 (.D) = J \Xs'>- (t) ( x ) - Xs(t)(x))dx, By Lemma [33] and by 
Theorem 12.11 the signed distance functions converge in C 2 ' a (0) uniformly in time and therefore 
this term converges to zero. The convergence of the term involving Xs^-(t) can be seen in the 
same way. □ 
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